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Abstract 

The Thomas- Fermi screening of non-Abelian gauge fields by fermions or screening of 
gluon fields in quark matter is discussed. It is described by an effective mass term which 
is, as with hard thermal loops, related to the eikonal for a Chern-Simons theory and the 
Wess-Zumino-Witten action. 
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The screening of electromagnetic interactions in a plasma of charged particles, the 
so-called Debye screening, is quite well-known. An analogous effect, the Thomas-Fermi 
screening, occurs in the case of a degenerate charged fermion system such as the electron 
gas in metal [1]. Both effects can be easily understood, in a field-theoretic language, by 
calculating the one loop photon propagator in which the charged particle propagators are at 
finite temperature and density. (The two effects are in fact quite similar; for the degenerate 
fermion gas, the excitations generated by the propagating photon are particle-hole pairs 
which behave like plasma background). The non-Abelian analogues of these screening 
effects are of considerable interest especially in view of the possibility of obtaining hot 
and dense quark matter systems in heavy ion collisions [2]. A non-Abelian Thomas- Fermi 
effect can also be of interest in calculations of the equation of state for quark matter inside 
neutron stars. Non-Abelian Debye screening and related effects have been intensively 
investigated over the last few years [3-8] . It is well-known that a proper calculation of this 
effect involves summing up all the hard thermal loop Feynman diagrams [3] . This results 
in a gauge-invariant, nonlocal, effective mass term for the gauge bosons (or gluons in a 
chromodynamic context). This mass term has many nice properties being closely related 
to Chern-Simons and Wess-Zumino-Witten (WZW) theories [5,6]. 

In this letter we analyze the non-Abelian Thomas-Fermi screening for degenerate 
quark matter of finite baryon number. Since the quark contribution to the two-point 
function for gluons is similar to that for photons, there should be Thomas-Fermi screening 
for quark matter. For reasons of non-Abelian gauge invariance, as with Debye screening, 
there will be higher point contributions, the whole series again summing up to an effective 
mass term. This term will have the same structure as the effective action for hard thermal 
loops; the numerical value of the screening mass, however, will be determined by the 
chemical potential rather than temperature. 

Consideration of the one loop two-point function shows that the screening mass ~ g[i 
where g is the coupling constant and \i is the chemical potential. It is then clear that 
a higher loop diagram in which such a term is inserted can give contributions of the 
same order for the integration range of loop momenta ~<7/U. This is exactly as in the 
hard thermal loop case. One must therefore sum up diagrams with loop momenta ~ 
and external momenta ~<7/U. This effective action must then be used for a self-consistent 
evaluation of the screeening mass. We obtain this effective action in what follows. 

Let us start by considering the two-point function for gluons. The conserved baryon 
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charge is given in terms of the quark field q(x) by J q^q. With a chemical potential term 
fi J q^q added to the action, the quark propagator is given by 

S(x, y) = (Tq(x)q(y)) = [ -^^[{6(x° - y G )a pl • pe'^^ + p7 • P 'e^'^}+ 

J {2ir) Z P 

-e(y° - x°){/3 p7 • pe-^-v) + a vl ■ p'e ip '^-^}] 

(1) 

where p° = \p\, p = (p°,p), p' = (p°, —p) and 9{x) is the step function. Also 

a p = 1 - n p , (3 P = n p (2) 

The distribution functions n p ,n p corresponding to quarks and antiquarks respectively are 
given by 

- 1 _ _ 1 

U P ~ e (pO- M )/ T + P n P ~ e ( P "+M)/T + ! 

The one-loop quark graphs are given by the effective action 

T = -zTr log(l + S'j ■ A) (4) 

= —it a gA a is the gluon vector potential, t a are hermitian matrices corresponding to 
the generators of the Lie algebra in the quark representation. In the above expression for 
T a functional trace is implied as well as the trace over the spin and color labels. The 
two-gluon term in V is given by 

T^= l - J d*x d 4 y Tr[ 7 • A(x)S(x,y) 1 ■ A(y)S(y,x)] (5) 

Using equation (1) and carrying out the time-integrations we get 
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where 
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d 4 k d 4 k' (7) 



d/i(k) = (2n) 4 5^(k + k') 
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T(p, q) = Tr [7 • A(k) 7 • p 7 • A(fc') 7 • g] 
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and p = q + k in equation(6). 

The ie's can be taken to go to zero at this stage. They were introduced for conver- 
gence of time-integrations and contribute to the imaginary part. Here we are interested 
in screening effects which are described by the real part of the two-point function. (Also 
for many physical situations, the relevant imaginary part is that of the retarded function 
which is not directly given by the above time ordered function [6].) Further we are in- 
terested in a degenerate gas of quarks; it is therefore appropriate to consider T ji. As 
T/fj, — > 0, the antiquark occupation numbers n p — > (for positive y). Equation (6) then 
simplifies as 



r<2) = -U d '^^w [TM w^ +r(M V 
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As explained in the introduction, the relevant kinematic regime is \p\, \ q\ ^> \k\, so that 
p° -q° -k° w -jfe • Q, p° + q°±k ^ 2q°, Q = (1, q/q°). Further, 
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T(p, q) w 8g° • QA 2 ■ Q) 

Tip', q) w T(p, q') w 4g° • Q'A 2 ■ Q + A x ■ QA 2 ■ Q' - 2A X ■ A 2 ) 

where A\ = A(k), A 2 = A(k'), Q' = (1, —q/q )- Expression (8) now simplifies as 

r^-i/w^^^^^-t-^,^.^ (1Q) 

A x ■ QA 2 ■ Q' - 2A X ■ A 2 )] 

Using J d 3 q jjfpfiQ) = — J d 3 q and some properties of Q-integration, T^ 2 ) sim- 

plifies to 

r<2) =~\S / (S 8 tr ^ 2 - " if- 4 - 4 - " f| - 41 - 42 -] 
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where A + = A_ = and we have taken the limit of T — > (small compared to 
((/). As with hard thermal loops, we can transform this back to coordinate space and write 

r< 2 > = y rfO tf( 2 >(A+, A_) (12) 



where 



K(A + , A -)=[J d4x tr(A + A_) + mI(A + ) + mI{A_)] 
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J dQ defines integration over the orientations of Q. K^> in equation (12) denotes the 
terms in K which are quadratic in A's; z and z are the Wick-rotated versions of x ■ Q' and 
x ■ Q respectively and x T is transverse to Q, i.e Q ■ x T = 0. /(A_) is obtained from I(A + ) 
by z <-> z. I(A + ), as explained elsewhere, is essentially the eikonal function for Chern- 
Simons theory. K{A + , A J) can be related to the WZW-action for a hermitian matrix 
M^M defined in terms of A±. 

Consider now the three-point function. One has twenty-four terms with denominators 
involving different combinations of loop momenta and external momenta. The simplifi- 
cation of this expression proceeds in much the same way as for - one can neglect 
terms with denominators that are of the order of fx and leave the terms that involve dif- 
ferences between loop momenta, which are of the order of k <C fx. After that, the relevant 
contributions will be 

d 3 q r rip 1 n r 1 



J (2tt)3 1 



r (3) = / 2 I "p A. + 



k°+p°- q° r°-q -k° 3 q° - r° + k° 3 p° - r° - k° ^ 

where p, q, r are the loop momenta, fci, ki and k% are the momenta of external gluons and 

— * — * , . 

p = q — k2,r = q + ks. After performing the dq -integral the final expression is 

r<s> - ifK" 2 ^^ + ^ ■ ■ ■ « p»> 

Using definition (13) one can show that 

r( 3 ) = _£1 j dn kW(a+,a_) (16) 

The fact that the same coefficient appears in both (12) and (14) is crucial; this guarantees 
gauge invariance of the full effective action V. The non-Abelian gauge-invariant completion 
of + is given by the full K of equation (13). The final answer is thus 

r = -£^/ dQK(A +1 A.) (17) 



We now turn to the question of the coefficient of K in equation (17). In calculating the 
screening effects using T, we should not change the Lagrangian for the theory, which may 
be, say, chromodynamics; we should only rearrange terms in the perturbative expansion 
and sum certain classes of diagrams. There are then two different but closely related ways 
of proceeding. We write the action as 

S = S + A / dn K 

J , (18) 
So = S Q cd -m 2 dQK 

Sq is used to define propagators and vertices and A is treated as being nominally one 
loop order higher than Sq. In the first case, we use the value as calculated above ( or the 
analogous value for hard thermal loops) for to 2 . After calculating higher order corrections, 
A is set to m 2 ( so that S — > Sqcd )• The corrections are in general non- vanishing and this 
is useful if the corrections are small compared to the lowest order value. The alternative 
is to keep m 2 as an arbitrary parameter and choose A ( as a function of to, say A (to) 
) so as to cancel out all the corrections. Upon setting A to to 2 , we get a gap equation, 
A (to) = to 2 , which can be solved for to. In the case of the hard thermal loops, the first 
approach is satisfactory. The relevant distribution for the momentum-integration is 

/(»>*= < 19 > 

The lowest order calculation only evaluates the contribution from the region of g~T. The 
probability contained in the region q > \T for the distribution is approximately 0.99. We 
can therefore expect that the neglect of the low (/-regime is not very significant for the 
numerical value of the screening mass ( within a calculational scheme with resummations 
as explained above ). For the case of Thomas- Fermi screening, the relevant distribution is 

f(q)dq » Oiji - q)q 2 dq (20) 

The probability contained in the region q > ^ is now 0.875. We thus expect that the lowest 
order value of the coefficient of — f <ifi K in equation (15), viz., / u 2 /47r 3 , is somewhat less 
accurate than the analogous quantity for hard thermal loops. In this case a self-consistent 
calculational scheme for including higher order effects might be more appropriate. 

After this work was completed, we became aware of a paper by Cristina Manuel where 
non-Abelian Thomas- Fermi screening is discussed in a kinetic theory framework and results 
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similar to ours are obtained [9] . Our approach, based on evaluation of Feynman diagrams, 
is complementary to this work. We thank Cristina Manuel for bringing this work to our 
attention. 
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